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Abstract. We present a new method for determining the frequencies and wavef unctions 
of black hole quasinormal modes (QNMs) and Regge poles. The key idea is a novel ansatz 
for the wavefunction, which relates the high-Z wavefunctions to null geodesies which start 
at infinity and end in perpetual orbit on the photon sphere. Our ansatz leads naturally to 
the expansion of QNMs in inverse powers ofL = / + l/2 (in 4D) , and to the expansion of 
Regge poles in inverse powers of u>. The expansions can be taken to high orders. We begin by 
applying the method to the Schwarzschild spacetime, and validate our results against existing 
numerical and WKB methods. Next, we generalise the method to treat static spherically- 
symmetric spacetimcs of arbitrary spatial dimension. We confirm that, at lowest order, the 
real and imaginary components of the QNM frequency are related to the orbital frequency 
and the Lyapunov exponent for geodesies at the unstable orbit. We apply the method to 
five spacetimes of current interest, and conclude with a discussion of the advantages and 
limitations of the new approach, and its practical applications. 
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1. Introduction 

Quasinormal modes (QNMs) are damped resonances which play a key role in black hole 
dynamics [H El El H]. It was noted four decades ago [5] that a perturbed black hole 
will, after an initial response, emit radiation with well-defined frequencies and rates of 
damping, in the manner of a bell sounding its last dying pure notes [5]. Physically, a black 
hole quasinormal mode is a decaying resonance which satisfies a pair of causally-motivated 
boundary conditions (typically being ingoing at the event horizon, and outgoing at spatial 
infinity). Mathematically, 'quasinormal ringing' emerges from a sum of residues of poles 
of the Green function in the complex frequency domain [7]. Each pole corresponds to a 
quasinormal mode of a single complex frequency. The real part of the frequency corresponds 
to the oscillation rate and the (negative) imaginary part corresponds to the damping rate. 

In a spherically symmetric spacetime (such as the Schwarzschild spacetime), the QNM 
frequencies ui n are labelled by two integers: multipole / and overtone n > 0. The properties 
of the QNM spectrum depends only on the properties of the field (e.g. spin s) and the black 
hole geometry (e.g. mass M, charge Q, angular momentum J). 

Even after years of study, black hole QNMs retain a certain mystique. The frequency 
spectrum has a beautiful and ornate structure (see e.g. Fig. 5 in [2J). The asymptotic 
properties of the highly-damped frequencies of the Schwarzschild hole have inspired 
speculation that QNMs may be linked to, e.g., Hawking radiation [S], black hole entropy 
[9], and loop quantum gravity \±0\ (however, subsequent study of charged, rotating and 
higher-dimensional black holes has shown that any such links cannot be as simple as originally 
envisaged [HH2])- QNMs are also of relevance to holographic principles such as the AdS/CFT 
conjecture [TBI IT3] . Here we focus merely on the role that QNMs play in wave propagation 
on black hole spacetimes, in the weakly-damped regime I > 2, I > n. 

It has been known for many years that quasinormal modes are intimately linked to the 
existence and properties of unstable null orbits [TS1 H3 US]. In the eikonal regime (I ^> n), 
the real part of the complex QNM frequency is related to the angular frequency of the orbit 
Q. It was recently demonstrated [19] that, in any static spherically symmetric asymptotically 
flat spacetime, to leading order in I the QNM frequency is 

lim Re(uji n ) = Ql, lim Im(^ n ) = — {n + 1/2)|A|, (1) 

l^oo I— >oo 

where A is the so-called Lyapunov exponent: the (inverse of the) instability timescale of 
the unstable orbit [201 121] ■ Motivated by this idea, we show here that the QNM frequency 
for spherically-symmetric black holes may be written as an expansion in inverse powers of 
L = I + 1/2, as Uin — zu^L + Wq ^ + w\ ^ L 1 + w\ ^ L 2 + . . . We develop a simple method 
for determining the expansion coefficients {w^} which may be taken to very high order. A 
key advantage of the method is that it provides expansions (in powers of L^ 1 ) for the radial 
wavefunctions, which are well-defined over all radii r > r^. It is straightforward to adapt the 
method to find an expansion of the Regge poles [221 123] in powers of u, as we will show. 

This paper is structured as follows. In Sec. [2] we briefly review the existing methods for 
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computing QNM spectra. In Sec. [3] we introduce the expansion- in- L method by application to 
the Schwarzschild spacetime, and validate against existing numerical and asymptotic results. 
In Sec. [4] we adapt the method to find Regge poles of the Schwarzschild spacetime. In Sec. [5]we 
generalize the QNM analysis to static spherically-symmetric spacetimes of arbitrary spatial 
dimension, and confirm the result ([!]). In Sec. [6] we apply our method to five spacetimes 
of physical interest (Reissner-Nordstrom; Schwarzschild-de Sitter; Nariai; Schwarzschild- 
Tangherlini; the canonical acoustic hole). We conclude in Sec. [7] with a discussion of interesting 
applications of the method and work for the future. 

2. Existing Methods for Determining QNMs 

A wide range of numerical methods have been developed for determining QNMs. The following 
non-exhaustive classification was outlined in |24j: (i) Time domain methods, whereby the 
least-damped frequencies are extracted from long, stable time evolutions [SJ [25]; (ii) Direct 
integration in the frequency domain as pioneered in [26]; (iii) Inverse potential methods 
whereby the 'exact' potential is replaced with a simple potential whose spectrum is known [2*7] ; 

(iv) WKB methods adapted from standard techniques in quantum mechanics [2H1 EHJ [30], [31] ; 

(v) Phase-integral methods using integration along anti-Stokes lines in the complex-r plane 
[321 1331 ITT] ; (vi) Continued fraction methods, first employed in the 1930s to find electronic 
spectrum of the H + ion, and adapted to black holes by Leaver [31]. The relative advantages 
and disadvantages of these methods are discussed in several review articles [JJ [2J H]. 

For numerical calculations, the continued fraction method (vi) is fast, reliable and 
accurate, and with certain modifications |35J, is robust at large overtones n and angular 
momenta I. However, the method itself provides relatively little physical insight. In contrast, 
WKB methods [28l [29l [301 [31] provide insight at the expense of accuracy and applicability. For 
example, the Bohr-Sommerfeld rule 'explains' the link between the low-overtone frequencies 
and the shape of the peak of the radial potential [281 EH]- For improved accuracy, the WKB 
method may be extended to sixth order [31] and above; nevertheless it still breaks down at 
large overtone n ^> / or if the potential is complex. 

In this work we develop a new method which bears some passing resemblance to the 
WKB approach [281 ESI EQj EJJ . However, instead of matching together independent solutions 
across a (real) potential barrier [36J, we instead seek a single solution that is valid everywhere 
outside the horizon. It turns out that, once the correct geometrically- motivated ansatz is 
employed, the wavefunctions and frequencies may be expressed as expansions in the angular 
momentum parameter 

L = 1 + 1/2. (2) 
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3. The Expansion of Schwarzschild QNMs 

Without further ado, let us introduce the method by applying it to the Schwarzschild 
spacetime, perturbations of which are governed by the 'master' equation 



dr: 



f(r) 



(L 2 - 1/4) 2(3 



(3) 



where f(r) = 1 — 2/r and r* = r + 21n(r/2 — 1) and (3 = 1 — s 2 with \s\ = 0, 1, 2 for scalar, 
electromagnetic and gravitational perturbations [6], respectively. Note that we have fixed the 
black hole mass to M = 1. Solutions which are 'ingoing' at the horizon (as r* — > — oo) satisfy 



U Z( r ) ~ S .(in) 



+ A 



(out) 



r* 



-oo, 
-oo. 



(4) 



Quasinormal modes correspond to the (complex) frequencies uoi n at which A[™ = 0. 



3.1. Method 

Let us introduce the following ansatz 

Ulu(r) = exp \ iuo / 11 + 



1/2 



(5) 



Note that the integrand in the exponent goes as 1 + 0(r~ 2 ) as r* — > oo, and as — 1 + Oir — 2) 
as r* — > — oo. Hence QNM boundary conditions are automatically satisfied if Vi u (r) is well- 
behaved in both limits. The geometric motivation behind this ansatz is explored in Sec. [5j 
Substitution of ansatz (Jsl) into ^ leads immediately to 



d f ,.dv 
dr V dr 



+ 



2iu[l + 



6 



27u 2 - L 2 27iuj 
s + — ^- 



1/2 



-1/2 



dr 



1 

4^2" 



2/3 



(6) 



At first glance, Q looks more complicated than ([3]), but it has the advantage of being 
amenable to an expansion in inverse powers of L. To look for the 'fundamental' (least- 
damped) n = mode, let us try the expansion 

exp (of (r) + L-^f^r) + L~ 2 Sf (r) + iT 3 Sf (r) + . . .) (8) 

Upon substituting into ([6]), and collecting together like powers of L, we obtain a set of 
independent equations which may be used to determine the coefficients vJk and the radial 
functions Sjfe(r) (briefly dropping the superscript for brevity). We find 



ViJr) 



L 2 : 
L 1 : 



27zu 2 _ x -1 = 
2iw-x [ 1 + 



1/2 



= +1/V27 

^0 + ^ + 5"^ ( 1 + 



-1/2 
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+ 







1/2 



1 + 



r 

-1/2 



J (ro-iS'J 



S' o ) + f{£% + (S>) 2 ) + f'S' + 



1 



etc., where ' denotes differentiation with respect to r. The QNM frequencies are found by 
imposing a continuity condition on {S' k (r)} at the unstable null orbit radius r = 3. That is, 

3 which implies that 



we demand that S' is continuous and differentiate at r 



w, 



(0) 



—I 



2V27 



and 



dS, 



(0) 



27 



dr 2(r + 6)(r-3) 
We may continue in this way, iteratively 



1 + 



6 



1/2 



(13) 



(14) 



First, we write down the equation for the 



(1 — k)th power of L, which contains the unknowns w k and S^'' . Next we insist that S' k 
is regular at r = 3, to fix tu^. Then we rearrange to find S^' , and differentiate once to 
obtain S^" which features at the next order. If we wish to construct the wavefunction we 
also integrate sj®' to obtain S^\r) (whereas if we only require the frequency expansion this 
step may be skipped). Finally we move to the next order, sending k — > k + 1, and repeat 
the steps. It is straightforward to implement the procedure in a symbolic algebra package to 
extend the expansion to arbitrary order. The lowest-order coefficients for perturbations of 
spin \s\ = 0, 1, 2 of the Schwarzschild spacetime are given in Table [TJ The expansion may be 
carried to 15th order and beyond with relative ease. For instance, the frequency of the n = 
gravitational (s = —2) QNM is 

V27u l0 — L — 0.5i - 1.30092593L- 1 + 0.2046039HL- 2 - 0.56376775L~ 3 

+ 0.254543922L" 4 - 0.19978348L~ 5 - 0.05688153iL^ 6 + 0.14611217L^ 7 

- 1.16950214iL~ 8 + 0.04477903L" 9 - 3.49565902iL^ 10 - 1.50293894L^ n 

- 6.56185182iL~ 12 - 4.40093979L~ 13 + C(L" 14 ) (15) 

(expressing rational coefficients to 8 decimal places). We note the following: (i) the method 
provides a neat expansion in inverse powers of L = 1 + 1/2, and a globally- valid (in r) 
expansion ^ of the QNM function; (ii) odd (even) powers of L in the Schwarzschild QNM 
frequency expansion ([7]) have coefficients which are purely real (imaginary); (iii) the series is 
probably not convergent for low L < 1, rather one expects that Eq. ((151 is an asymptotic 



7(0)' 



expansion. That is, for a given L, series (15) may be formally divergent; nevertheless excellent 
approximations to the frequencies may be obtained by truncating at finite order k m such that 



\zu 



(0) 



< 



(0) 



The method can be extended higher overtones with the ansatz 



VinVr) 



i=l j=l 



1 - 



exp 



L~ 1 S\ I 



(»), 



r + 



(16) 
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Table 1. Series expansion coefficients Wk = \/27m k for the fundamental (n = 0) 
Schw. QNM frequencies, in expansion -v/27 Muo^ n=0 — Y^k=-i w kL~ k , for scalar (/? = 1), 
electromagnetic ((3 = 0) and gravitational ((3 = —3) modes. 





L 1 




L- 1 


L- 2 


L" 3 


L- 4 


L~ 5 


L~ 6 


n = 1 
n = 2 
n = 3 


1 
1 
1 


-H 
- 7 -i 


341 


1121 ■ 
2592 
905 • 
7776 
8603 • 


538673 


63615613 ■ 


23221986259 


453971786399 • 


216 
461 


629856 
146381 


120932352 
115304885 • 


39182082048 
19631537719 


940369969152 
10658412628805 ■ 


216 
641 


157464 
289949 


362797056 
31522981 • 


39182082048 
125344410311 


8463329722368 
49100100676571 • 


216 


7776 


1259712 


362797056 
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Table 2. Series expansion coefficients — y/Tfzu}! 1 for gravitational (s = —2) Schw. QNMs 
of higher overtones (n = 1,2,3) in expansion \/TJu)i n = w kL~ k . 



At the zth step of the interative procedure we require continuity of the first n — 1 non- 

(n) 

vanishing derivatives at r = 3 to determine the a\j , and nth to determine the correction to 
the frequency, tu t - n \ When these conditions are imposed, we are left with an explicit equation 
for S; n ^(r). 



3.2. Results: QNM Frequency Expansion 

Expansion coefficients for the frequency of the 'fundamental' mode are listed in Table [TJ 
for perturbations of spin (scalar), 1 (electromagnetic) and 2 (gravitational). Expansion 
coefficients for the n = 1,2,3 quasinormal mode frequencies of the gravitational field on the 
Schwarzschild spacetime are given in Table [2] 

The lowest expansion coefficients for arbitrary spin (3 = 1 — s 2 and arbitrary overtone 
number n can be written as 



27zu 



(n 



27zu 



(n 



27vu\ 



(n 



21 w 



21 w 



21 w 



(n 



(n 



(n 



1 

-iN 
P 5N 2 
3 ~ l56~ 



-iN 

-?- + 

21 
21 



9 



+ 



115 
432 
235iV 2 

3888 



1415 
15552 



204iV 2 + 211 „ 854160iV 4 - 1664760A^ 2 - 776939 



iN 



+ 



3888 
1100 iV 2 



46656 



+ 40310784 

2719 „ 11273136AT 4 - 52753800A^ 2 + 66480535 

13 + 



2902376448 



(17) 
(18) 

(19) 
(20) 

(21) 
(22) 



where iV = n + 1/2. 
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I = 2,n = 


I = 2, n = 1 


I = 3,n = 


I = 3, n = 1 


Ctd. Frac. 


0.373672 -i0.088962 


0.346711 -i0.273915 


0.599444 -i0.092703 


0.582644 -i0.281298 


12th order 


0.373679 -i0.088934 


0.346831 -i0.273844 


0.599443 -i0.092703 


0.582645 -i0.281297 


6th order 


0.373642 -i0.088671 


0.347895 -i0.272766 


0.599439 -i0.092684 


0.582713 -i0.281206 


WKB (6th) 


0.3736 -i0.0890 


0.3463 -i0.2735 


0.5994 -i0.0927 


0.5826 -zO.2813 



Table 3. Numerical Accuracy of the Gravitational QNM frequencies. Comparing the results 
of the new expansion method at 6th and 12th order (with truncation errors of order 0(L~ 6 ) 
and (D(L~ 12 ), respectively) with numerical results from the continued fraction method [32] 
and 6th order WKB results [31]. 



3. 3. Validation 

In Table [3] the results of the expansion in powers of L~ x are compared against the continued- 
fraction results of [31], and the WKB results of [3T], for Schwarzschild gravitational modes 
with / = 2, 3 and n = 0, 1. At 6th order (i.e. including all terms in series Q up to and 
including tUj L" 5 , resulting in truncation error 0(L~ 6 )), the accuracy of the new method is 
comparable with the WKB method. It is simple to extend our method to, e.g., 12th order to 
obtain greater accuracy, as we show in Table [3] 



We checked the series expansion coefficients (17-22) for general spin j3 = 1 — s 2 and 
overtone n against the WKB results given in Eq. (3.1) of |30j. We find agreement up to and 
including coefficient m% (power L~ 2 ). This is expected, because the truncation error in the 
WKB results of [30] is of order 0(L~ 3 ). 

The accuracy of the expansion is examined in Fig. [T] and [2j Figure [T] compares the QNM 
frequencies found via the expansion method with the results of the continued fraction method 
[34] . As expected, we find good agreement if I > n, but not always in the regime I < n. In 
general, the expansions of the higher-n modes are less accurate. This is clearly demonstrated 
in Fig. [2J which shows the difference between the series expansion (taken to orders L~ 5 and 
L -10 ) and the exact frequencies, for the modes n — . . . 5. 

Figure [3] shows the magnitude of the coefficients in the series expansion of the frequencies. 
It shows that the coefficients do not seem to tend to zero as k — > 0, though they do become 
small at certain intermediate orders. This supports the view that Eq. Q is an asymptotic 
expansion. 

3.4- Results: Outgoing coefficient and the wavef unction 

A key advantage that the expansion method holds over many other approaches (see Sec. [2]) 
is that it furnishes us with simple approximations for the wavefunctions [Eq. ^ and (16)] 



which are valid at all r > 2. Here we present the results for the outgoing coefficient in 
a compact form, and examine the shape of the wavefunction. 

The integral appearing in ansatz ^ can be evaluated without much difficulty. Let us 
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Exact Schwarzschild QNM frequencies (+) compared with l_arge-(l+1/2) Expansion (x and lines) 
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Figure 1. Gravitational QNM frequencies of Schwarzschild black hole. The exact QNM 
frequencies (+ symbol, determined via the continued- fraction method [34]) are compared 
with the results of the expansion in L^ 1 (x symbols, dotted lines) taken to order (D(L~ W ). 



define 



X(r) = exp ( iu / [(1 + 6/r) 1/2 (l - 3/r)] f^dr 



(2 + V3Y 



Siu) 



1 + X 

2-x 



exp l iuj 



r(l — x) 



where 



x 



1 + 



1/2 



The outgoing coefficient defined by Eq. (4) may be computed via 



.4 



(out) 
In 



("2) ? 



.i + ELiEr=i(-2) fe 4?^ 



a 



>ut / (n) T -k 

^exp 2^71 'L 



vfc=0 



where we have defined the quantities 

C m = lim ^(r)/e~ iwr * = exp (iw[6 



r->2 



27 + 81n2-61n(2 + 



C out ee lim X(r)/e +iu)r * = exp (zu;[3 - Vff + 4 In 2 - 6 ln(2 + V3)} 

r^oo \ 



lim S^\r) 

r— >oo 



and we choose the constant of integration such that lim r ^2 S 

7 (n) 



(«)/ 



(23) 
(24) 

(25) 

(26) 

(27) 
(28) 
(29) 



= 0. The leading-order 

function S^'fr) is independent of the spin of the perturbing field, and may be written 



On an Expansion Method for Black Hole Quasinormal Modes and Regge Poles 



9 



Numerical Error at 1 1th order [trunc err. 0(L' 11 )] 
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Figure 2. Accuracy of Expansion of Schwarzs child gravitational QNM frequencies in inverse 
powers of L = 1+1/2. The four plots show the difference between the exact QNM frequencies 
(continued-fraction method) and the approximations from the asymptotic expansion method. 
The frequency expansion lui u = X)fc= 

truncated at orders k m = 5 (left) and 
k m = 10 (right). The upper (lower) plots show the accuracy of the real (imaginary) part of 
the frequency on a logarithmic scale, for overtones n = . . . 5 (dotted lines). 



compactly as 



S ( o l \r) = ^ ln(2/x) + 2{n + 1/2) In 



2 + v^' 



(30) 



Note 5*0 (r) is purely real, hence it tells us about the amplitude of the QNM in the large-L 
limit. The leading-order coefficient defined in (29) is therefore 



(n) 

7o 



1/2) In 2 + (n + 1/2) In ([2 + V?>\/2\ . 
Hence, at leading order in L, the outgoing coefficient is simply 



a (° ut ) 



-2) -Q- exp \^f K 



n+l/2 



-icj irl (3+41n2) 



(31) 
(32) 



for all spins. 



Calculation of (26) is straightforward for the 'fundamental' mode, for which we find 

(2 + Vz) 1,2 e- 1 ^ ( 3 + 41n2 ) e xp ( f^lfL-" ) (33) 



A ( out ) 
yl Z,n=0 



,fc=l 
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Figure 3. Coefficients in the Expansion of the Fundamental Schwarzs child QNM Frequency 
(Scalar Field). The plot shows the coefficients vj 1 -- ^ in the expansion = 2fcL-i ' '" 1 



L~ 

The lower plot shows the same data as the upper plot, using a logarithmic scale. [N.B. The 
coefficients are known as rationals (i.e. exact)]. 



where the first few coefficients for the fundamental mode are 
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Higher-order approximations can be found by computing the functions S 



in). 



(34) 
(35) 
(36) 
(37) 
and the 



coefficients . Typical radial profiles of the functions Sj^\r) are shown in Fig. ^ We 
see that the functions start from zero at r = 2 and tend towards a constant as r — > oo, 
without any significant oscillation between these limits. 

Figure [5] shows examples of the QNM wavefunctions. Since the frequency has a negative 
imaginary part, it is no surprise to find that the wavefunctions diverge as r — > 2 and as r — > oo. 
This makes it rather hard to gain physical insight from examining the radial functions alone. 
To assess the QNM contribution to physically-realistic scenarios, we need to compute the 
'QNM excitation coefficient' [13]. This quantity is obtained from the residues of the poles 
in the Green function at QNM frequency. In a forthcoming work [37] we will show how to 
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Radial Functions in Wavefunction Expansion v(r) = exp(S (r) + L" 1 S-|(r) + ... ) 
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Figure 4. Radial profile of the functions S k (r). Note that functions with k even (odd) are 
purely real (imaginary). The functions for higher overtone n are similar in nature. 




Figure 5. Gravitational QNM wavefunctions of lowest overtones (n = 0, 1) at L = 1 + 1/2 — 
4. Note that the wavefunctions diverge exponentially in both limits, — > ±00. 



, by 

UJ=U} ln 



calculate the residues of the poles, and in particular the key quantity dA^ /doo 

combining the expansion method with standard WKB techniques. 

An alternative set of physically-relevant wavefunctions that do not diverge in the limits 
r — > 2 and r — » 00 are the so-called Regge poles (shown in Fig. [6J, to which we now turn our 
attention. 

4. The Expansion of Schwarzschild Regge Poles 

For a given angular momentum I, there are an infinite set of QNMs with frequencies uj\ n 
defined by A^y = 0. Conversely, for a given frequency uj there are an infinite set of Regge 
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poles, corresponding to complex angular momenta X un = + 1/2, defined by the condition 
A w = 0- Regge poles are a key concept in the so-called 'Complex Angular Momentum' 
(CAM) method |38j, which has been applied in a variety of contexts. For example, Regge 
poles received much attention in high-energy physics in the 1960s [39J. In the context of black 
hole physics, it has been shown that the scattering amplitude for a monochromatic planar 
wave impinging upon a black hole can be split into a background integral and a sum over 
Regge poles 001 HI]- The background integral (approximately) corresponds to the 'classical' 
scattering cross section, and the Regge poles correspond to regular oscillations in intensity-vs- 
angle plot caused by the interference of parts of the wavefront travelling close to rays passing 
in opposite senses around the black hole (see [32] for a short review). 

Many of the methods used for finding QNMs can also be used to find Regge poles [221 23] . 
It is perhaps no surprise to find it is straightforward to expand Regge poles X un in inverse 
powers of to. We simply insert the expansion 

A wn = \™u + A< n) + X^lu- 1 + X^lu- 2 + ... (38) 



(r) 



( X " I) n + E E b ^~ j ( X " I) ^ ex P + ^T[ n \r) + ...) (39) 



into the differential equation ([6]), and group like powers of to, to obtain a set of equations. 



We then solve order-by-order to determine the coefficients X^\ bffi and radial functions 



Ti. , as before. We find that the first few expansion coefficients are 



k 

» 



\™ = V27 (40) 
A< n) = iN (41) 
a(b) = 60iV 2 -144/?+115 
1 432^ 



_ /1220A?- 6912/3 + 5555\ 

2 V 419904 J K 1 
(n) _ 18/3 2 [2316iV 2 + 479]/3 2357520iV 4 + 19382280iV 2 + 2079661 

3 ~ 6561^ 104976^27 1088391168\/27 

Xf = iN 



8/3 2 5/3[3716iV 2 + 2291] t 144920784iV 4 + 1871793480iV 2 + 593617841" , ( _, 



19683 17006112 2115832430592 

where N = n + 1/2 and (3 = 1 — s 2 . We have verified that the expansion in inverse powers of 
u) given in Eq. (40 45) is fully consistent with the WKB result recently obtained by Decanini 
and Folacci ([23], equation (12) and (16)), up to and including order uj~ 2 (but not beyond). 
As in the QNM case (c.f. Sec. 3.3) this is consistent with the expected truncation error in 
the WKB results. We may take the expansion to higher orders if desired; for example, the 
'fundamental' n = Regge pole for the gravitational perturbation is 

^ r— i 281^ 6649i 1044601^ 926224193i 
W) „=o ~ V u + - + ig44cJ - 209952c; 2 ~ 153055008c; 3 + 264479053824c; 4 

184851431845^ 71361067332161z 14390928366797903\/3 



257073640316928c; 5 166583718925369344c; 6 161919374795459002368c; 7 
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Re(A aJjn=0 ) Im(Aa; in= o) 


Re (error) Im(error) 


Expected truncation 
error (approx.) 


Moj = 2.5 
Exact 
15th ord. 
8th ord. 
WKB 


13.0897823 0.4950209* 
13.089782327 0.495020871i 
13.089782322 0.495020835i 
13.0899530 0.4950172i 


(0.0%) (0.00002%) 
(0.0%) (0.00002%) 
(-0.0013%) (0.0007%) 


(2. x 10~ 12 + 3. x 10~ 13 i) 
(2. x 10~ 7 + 2. x 10" 6 i) 


Moj = 1.0 
Exact 
15th ord. 
8th ord. 
WKB 


5.4358037 0.4714533i 
5.4358037 0.4714533* 
5.4357860 0.4714045i 
5.4381254 0.4713563i 


(0.0%) (0.0%) 
(0.0003%) (0.01%) 
(-0.04%) (0.02%) 


(3. x 10~ 7 + 1. x 10~ 7 i) 
(1.5 x lO" 4 + 4.3 x 10" 4 i) 


Moj = 0.5 
Exact 
15th ord. 
8th ord. 
WKB 


3.0310943 0.4112060i 
3.0319568 0.4106489z 
3.0243829 0.4019403* 
3.0434848 0.4109477i 


(-0.028%) (0.135%) 
(0.22%) (2.25%) 
(-0.41%) (0.06%) 


(1.8 x 10- 4 + 2.7x 10~ 3 i) 
(2.0 x 10" 2 + 2.7 x 10- 2 i) 



Table 4. Numerical Validation of Regge Pole Expansion Method. The Regge pole expansion 



Eq. (38) is compared against the results of the continued-fraction method ('Exact') and the 
WKB method obtained by Decanini & Folacci [53] . The accuracy of the expansions of A WTl at 



orders i 



and 



are compared. Expected 'truncation errors' (found from the magnitude 



of the last terms in the series) are stated in the last column. 



+ 



7717840397981223883i 



+ 



14117571610293670714747\/3 



139898339823276578045952c; 8 1223830676774023504745988096c; 9 



+ . . . (46) 



Note that the expansion (46) only works in the regime Moj > 0.5, but is not appropriate or 
effective near oj ~ 0. 

In Table [4] the results of the expansion method are validated against the (approximate) 
WKB and the (exact) continued-fraction results presented in [23]. We find that the high-order 
expansions are in general more accurate than the WKB results, at large and intermediate oj. 



However, at small oj < 0.5 the expansion (38) does not converge well, hence the results are of 
limited use. The WKB method [23], by contrast, gives results which seem relatively robust 
down to oj ~ 0.25. 

Wavefunctions of the lowest overtones (n = 0, 1, 2) at oj = 3 are shown in Fig. [6] Note 
that, unlike the QNM wavefunctions (JSJ), they do not diverge in the limits r* — > ±oo. 



5. Expansion of QNMs in Spherically-Symmetric Spacetimes 



Having established the utility of the method on the Schwarzschild spacetime, let us now seek 
to generalise the analysis. In Sec. 5+] and 5^ we derive the leading-order terms for the QNM 
frequencies in a general static spherically-symmetric spacetime. The analysis highlights the 
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Figure 6. Regge pole wave functions of lowest overtones (n — 0,1, 2) of gravitational 
perturbations of a Schwarzschild black hole at Moj = 3.0. 

link between the geodesic equations and the expansion method. In Sec. [6] we look more 
carefully at some specific cases: Reissner-Nordstrom; Schwarzschild-deSitter; Nariai; higher 
dimensional black holes; and the canonical acoustic hole. 

5.1. Critical Orbits 

Let us assume a (d + 2)-dimensional static spherically-symmetric spacetime with line element 

ds 2 = -f(r)dt 2 + g~\r)dr 2 + h{r)dQ 2 d) (47) 

where dQ 2 * is the natural line element on S d and we assume that h(r) > everywhere. This 
includes all the spacetimes mentioned above, indeed for all of them except Nariai we may 
take h(r) = r 2 while for Nariai h(r) = 1. 

Null geodesies on this spacetime may be taken to lie in the equatorial plane where the 
orbital equation may be written as 
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Here the constant b = L/E may be interpreted as an 'impact parameter', where L = h(r)4> 
and E = f(r)t are constants of the motion and ' denotes differentiation with respect to an 
affine parameter. 

Let us define a new function 

i m 

b 2 h{r)' 

Now assume that there is a critical impact parameter b = b c , for which k 2 (r, b c ) has a repeated 
root so there exists some pair of constants {r c , b c } such that 

dk 2 



k 2 (r, b) 



(49) 



k (r c ,b c ) 







and 



dr 



(r c , b e ) = 0, 



(50) 



Assuming that the repeated root is a double root let us then write 

k c (r) = sgn(r - r c )\/k 2 (r, b c ) = (r - r c )K(r), (51) 

so k c (r) is positive for r > r c and negative for r < r c . 

It is conventional to study the stability properties of this orbit by writing the radial 
equation as 

• 2 L 2 g(r) 



f(r) 



-k 2 (r,b) = V r (r) 



Then the period of the circular orbit is given by 2ir /VL C where 

a 



7(r c 



7'(r c ) 



t 



h(r c ) V h '( r c) 



l/b c 



and the instability timescales are determined by the Lyapunov exponent [20 



A = = V9(rc)h(r c )K(r c ). 



(52) 



(53) 



(54) 



5.2. Wave equations 



To motivate our form of the wave equation we start by considering a scalar wave equation 
with potential (□ — V)<p = 0. Using the fact that the eigenvalues of the Laplacian on S d are 
-1(1 +d-l) (with multiplicity ( d f) - ( d+l ~ 2 )), the wave equation separates to give a radial 
equation of the form 



fgh 



-d/2 



d 

dr 



fgh d ' 2 



dR\ 



- 



dr J 



+ 



M 

h(r) 



(L 2 -(d-l) 2 /4)-/F(r) 



R = 



(55) 



where L — I + (d— l)/2. Introducing f * by y 1 = l/y/Jg and changing dependent variable to 



R — h d / 4 u, this may be rewritten as 



d 2 u 
dr 2 



+ 



m 

h(r) 



(L 2 — (<i — l) 2 /4) — fV e s(r) 



u 



where 



V eS (r) = V(r) + h- d / 4 f- 1 ^fg 



d 
dr 



dh d l A \ 
dr J 



(56) 



(57) 
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To provide a framework of sufficient generality we shall now study Eq. (56) but allow V e g(r) 
to contain L dependence at lower orders than L 2 , so 

Kfl(r) = V^(r)L + V (r) + V^L' 1 + .... (58) 

This form is motivated by, for example, gravitational perturbations of the Reissner-Nordstrom 
spacetime which we study below. 

Next we introduce the corollary to ansatz ([5]), namely, 



u(r) = exp I iu I b c k c (r)dr* I v (r) 



(59) 



where k c (r) was defined in Eq. (51). In all the cases we shall assume that we have a horizon 
as r* — > — oo, so fir) — > and correspondingly b c k c (r) — > —1. In addition as r* — ► oo we 
either have an asymptotically flat region or a cosmological horizon, so f(r)/h(r) — > and 
b c k c { r ) —> +1- Our ansatz thus encapsulates the desired ingoing and outgoing boundary 
conditions for the wave. Note that our assumptions at this stage have excluded, for example, 
Schwarzschild-AdS space-time. 



We now substitute ansatz (59) in to the radial equation (56) to obtain 



2ib c ujk c 1 r 



+ 



1 - b 2 c k 2 c )u 2 -t(L 2 -(d- l) 2 /4) - fV cS + iub t 



dk c 
dr 







(60) 



Once again, we expand to and v (r) in inverse powers of L, using Q and @. Making use of 
the key property 1 — b 2 k 2 (r) = b 2 f(r)/h(r) we balance at order L 2 for all r by taking 

' ' (61) 



a. 



Next at order L we require 

2iVfg(r - r c )KS' + 2(1 - b\{r - r c ) 2 K 2 )Q c u - fV^ + i^/fg [K - (r - r c )K'} = 0. (62) 



Setting r 



requiring S' to be regular there, we find 



fcV- 



lc 



2Q r 



«2 V 'QcKKc 



fcV- 



lc 



2a 



2 



(63) 



where a subscript c on a function of r denotes that function evaluated at r = r c . Note that, 



to leading order, the imaginary part of the frequency is set by the Lyapunov exponent ( 54 ) 



as anticipated in Eq. ([TJ (see also [HI [201 |2T] ) . Substituting back into Eq. (62), we may write 
an explicit and manifestly regular equation for 5*q: 

2iy/TgKS' Q = A[fV^} - iA[^TgK} - i^JgK' + 2(u /Q c )(r - r c )K 2 , (64) 

where we have introduced the notation 

F(r) - F(r c ) 



A[F](r) 



r c 



(65) 
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This process may readily be extended to higher orders yielding 

p 

^p+i = {-fc9c{Sp C + s icS' p - lc ) - i{u P /sic)^Jg c K c - \fJg c (\fJg)' c S'pc 



i=0 



- UjUp-i + fVpc - ^ ^ — —8 pi] )/(2il r 



i=0 



and 



2WfgKSL 



p+i 



-A 



fg(s; + S'iS'^)] - Aiy/fgiVfgys',] - i(u; p /n c )A[^/fgK} 



-A 



f(d-l) 
Ah 



i=0 



5 pQ + A[fV p ]-i(u p /n c )^fgK' 



p p 
2i sffgK ^ QJiS' v _J Q c + (2Q c u p+1 + J^w^/ft 

i=0 i=0 



, (V r - r c )K . 
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(66) 



(67) 



6. Examples 

To further demonstrate the utility of the above approach, we now apply to five cases of current 
interest. 



6.1. Reissner-Nordstrom Black Hole. 

The expansion method may be applied to deduce the QNM frequencies and wavefunctions 
of the charged (Reissner-Nordstrom) black hole, with a charge-to-mass ratio q. The relevant 
radial equation is P H [Ml IS] 



u 2 - f(r) 



L? - 1/4 «| S | q 2 r)\ s \ 



where 



+ 



Scalar 



u = 



K 



(68) 



(69) 



\s\ = \ 3 - a/9 + 4q 2 (L 2 - 9/4) |s| = 1 Electromagnetic 
3 + A/9 + 4g 2 (L 2 -9/4) |s| = 2 Gravitational 

and 77| s | = {2,4,4} for \s\ = {0,1,2} and f(r) = 1 — 2/r + q 2 /r 2 . For |s| = 1,2, we may 
expand k\ s \ in inverse powers of L, e.g. k\ 8 \ = ±2qL + 3 + 0(L~ r ). The resultant terms at 
order L in the radial equation mean that the odd (even) powers of L are no longer purely real 
(imaginary). 

Applying the method, we find the 'fundamental' (n = 0) gravitational mode is 

1/21 



L - 



2q i 
+ 



2a 



a + 3 2 V a + 3 



;i + «)(ll-25a + 84a 2 ) . 3V2q(a 2 - 1) ) + ^ ( _ 0) 



64(a + 3)a 2 L 



8(3a + a 2 ) 3 / 2 L 
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where b c = (3 + af/ 2 /{2(a + I)) 1 ' 2 and a = ^9 - 8q 2 and q is the charge-to-mass ratio of 
the black hole. Note that the formula breaks down, not when the twin horizons merge and 
disappear at q = 1, but rather when the unstable circular orbit disappears, at q = 3/\/8. 
More accurate numerical results for I = 2 are given in, e.g., [HI I4"5] . 
For the critically-charged case, q — 1, we obtain 

L 2 + \/2i 35 560 - 131\/2i 
4 16 256L 8192L 2 

+ 0(L- 5 ) (71) 



<^.n=0 



5031 - 2096 v^i 392512 - 137605 v^i , ^ , _,, 



131072L 3 16777216L 4 
The accuracy of the eigenvalues obtained from our series is competitive with the most 
accurate values in the literature [46]. For example, at 10th order we obtain uj^q = 
0.9657626 — 0.08700129/ which agrees with the I = 4 mode presented in [46] to the five 
significant figures given there. 

6.2. Schwarzschild-de Sitter Black Hole. 

The expansion method may also be applied to non-asymptotically flat spacetimes. For 



example, in Schwarzschild-de Sitter (SdS) spacetime, the line element takes the form (47) 
with f(r) = g(r) = 1 — 2/r — Ar 2 /3 and h(r) = r 2 , where < A < 1/9. The 'fundamental' 
(n = 0) QNMs of the scalar field in SdS are given by 



27 1 7-61?7 137 + 186877 - 2005r/ 2 

■ vi in = = L - 1- + 



1-9A l '" =u 2 216L 7776L 2 

5230 + 440043^ - 1274856?? 2 + 750851^ 3 
+ 2519424L 3 ^ ' 

. 590983 + 37791548// - 243504102// 2 + 340616636r? 3 - 135495065r/ 4 _ 5 
+% 362797056L 4 + ^ ' 

where r] = 9A. It is straightforward to extend to fields of higher spin, or to higher modes. 
The results compare well with, e.g., the 6th-order WKB results of Zhidenko ([Ej, Table 1). 

6.3. Nariai Space-time. 



Nariai space-time [H2] is the natural metric on the Lorentzian version of S 2 xS 2 , corresponding 
to /(r) = g{r) = 1 — r 2 and h(r) = 1. It provides a valuable check on our ansatz since 
the QNM frequencies can be determined exactly; for a scalar field with curvature coupling 
V = V eS = £R = 4£ 



Uln = ^L 2 + Ai-\-i(n+\) . (73) 

Our expansion precisely yields the large L expansion of this exact result, so at least in this 
case our expansion is not merely asymptotic. 
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6.4- Higher- Dimensional Black Hole. 

An uncharged spherically-symmetric higher-dimensional black hole in D = d + 2 dimensions 
(i.e. one dimension of time and d + 1 spatial dimensions) is known as a Schwarzschild- 
Tangherlini black hole jlH], and its metric is 

ds 2 = -f(r)dt 2 + f-\r)dr 2 + r 2 dQ 2 [d) (74) 

where f(r) = g(r) = 1 — (r/ l /r) a!_1 , and is the outer horizon. There is a circular orbit at 
r = r c with critical impact parameter b = b c where 

d + iV'w , , fd + i\ 1/2 , s 

— 2~ I r h and b c = I — — 1 r c . (75) 

The gravitational perturbations are divided into three classes, labelled 'tensor', 'vector' and 
'scalar' [49l [50] . The tensor perturbations are the easiest to analyse; they satisfy the radial 
equation 



j f(>) ( L 2 -(d-l) 2 /A + fd + d (d- 2)/(r) 



u = (76) 



r 2 ' 2r ' 4 

where recall L = Z + (d — l)/2. The lowest-order (in L) terms in the 'tensor' QNM expansion 
are 

+ (77) 

where we have used the result r 2 k' c {r c ) = b c \J d — 1. This result is, unsurprisingly, consistent 
with e.g. Eq. (51) in 1 1 9 and Eq. (12,13) in [31]. It should be possible to apply the method to 
obtain higher-order terms, and it may be possible to analyse 'vector' and 'scalar' perturbations 
in the same way, although the potential for the 'scalar' case is not positive definite, and it 
has more than one extrema [4~9l |50j [51] . 

6.5. Canonical Acoustic Hole. 

Black hole analogues - systems with some of the key properties of black holes - have received 
much attention in recent years [531 G3]- One of the simplest models for an analogue in 
fluid mechanics is the 'canonical acoustic hole', first proposed in [HS]- An 'horizon' forms 
where the bulk speed of the fluid (inwards) exceeds the speed of sound in the fluid. This 
raises the interesting possibility that QNM frequencies may one day be measured in the 



laboratory [5H]. The effective geometry of the 'canonical hole' is given by line element (47) 
with f{r) = g(r) = l — r^/r 4 and h{r) = r 2 (and d = 2). The 'fundamental' QNM frequencies 
are given by 

61 17 532843 4802843 rnfT -^ 

& c ov =0 = L-t- — -t - 25ig424L 3 + » 5668704^ + ° {L ] (?8) 

where b c = 3 3 / 4 /2 1 / 2 r/ l and L = 1 + 1/2, a result which is consistent with Eq. (53) in [SB]. The 



accuracy of the QNMs obtained via (78) is comparable or superior to the best estimates in 



the literature. For example, at I = 2, Eq. (78) gives uor h ta 1.4726 - 0.6087 ± [0.01 + O.Olz 
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which should be compared with ur^ m 1.41 — 0.70z found using a 6th-order WKB scheme 
[56] . A detailed study of the QNMs of acoustic holes is currently in progress [57] . 

7. Conclusion 

In this paper we have introduced a new method for obtaining QNM frequencies and 
wavefunctions through expansion in inverse powers of the angular momentum L (where 
L = I + 1/2 in four dimensions). We illustrated the method by applying to the Schwarzschild 
spacetime (Sec. [3]), to obtain an expansion of the QNM frequency (Eq. [7] and Eq. 17-22) 



that may be taken to very high orders (Eq. 15, Tables [T] and [2]) . We validated the frequency 
expansion against existing results (Table [3]), to show that it is accurate and rapidly convergent 
(Fig. [2J in the regime I > 2, I > n (but generally inaccurate/poorly convergent outside this 
regime, see Fig. [I]). In Sec. [4] we showed that the expansion method may also be applied to 
find Regge poles (Eq. 404451). In Sec. [5] we generalized the method to treat static spherically- 
symmetric spacetimes of arbitrary spatial dimension. The method was then applied to a 
selection of five such spacetimes of interest in Sec. |6j 

The expansion method complements the standard WKB approach [281 EHl EH [3TJ . In 
addition, we believe the method holds certain key advantages. It provides a simple high-order 
expansion for QNM (and Regge pole) wavefunctions that is convergent everywhere outside 
the horizon (in the regime of validity I > 2, 1 > n); it yields the complex coefficient 
(Eq. 32 and 33) which expresses the ratio of flux at infinity to flux at the horizon; and it 
provides additional insight into the link between orbiting null geodesies and QNMs. 

Let us briefly expand on this last point. The crucial step in the new method is a 



'geometric' ansatz for the wavefunction, given in Eq. ^ and Eq. (59). In Sec. 5.1 we showed 



how the ansatz is related to the family of 'critical' null geodesies that approach from infinity 
and end in perpetual orbit on the photon sphere at r = r c . In some sense, in this work we 



have conducted an expansion about these 'critical' geodesies. In Sec. |5.2| we proved a key 
result: in the eikonal limit / — > 00, the QNM frequencies depend only upon the properties 
of the 'critical' null geodesies, confirming the result of [19]. The real part of frequency is 
Re(a;; n ) = L/b c , where b c is the critical impact parameter (and b c = 1/Q C where fl c is the 
frequency of the null orbit at r = r c ), and the imaginary part is lm(uJi n ) = — (n + 1/2) |A|, 



where |A| is the Lyapunov exponent given in Eq. (54). Our analysis is valid in general for any 
perturbation of a static spherically-symmetric spacetime that satisfies 'radiative' conditions 
at both boundaries r* — > ±00 (a condition which excludes, for example, Schwarzschild-AdS). 

In the highly-damped regime n ^> l,n ^> 0, it is well known that QNMs have the 
asymptotic form ui n ~ 3? — inQ (see e.g. [63] for a discussion). Here 3? and 3 are real 
constants, where 3? depends on the geometry and the spin of the field, and 3 depends only 
on the geometry. For single-horizon asymptotically-flat black holes it has been shown that S 
is set by the surface gravity of the horizon [64]; for multi- horizon holes the situation is more 
complicated [65J. A possible challenge for the future would be to adapt the methods of this 



paper to investigate the highly-damped regime n>!,ti>0. 
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It is hoped that the existence and properties (e.g. mass, angular momentum) of 
astrophysical black holes may one day be inferred from observations of QNM ringing (in 
lightly-damped modes) at gravitational- wave detectors. We have seen that the relevant part of 
the QNM frequency spectrum is relatively easy to compute; however it is much more difficult 
to obtain an estimate of the degree of ringing excited by a given source. The theoretical 
framework has been in place for over twenty years [7] , yet because of the technical complexity 
only a few attempts at practical calculations have been made [7J |59j [601 ETJ [621 H3j. This 
remains an interesting avenue for further study (and the Regge pole formalism, developed in 
pEOl l38~l |2"2| [23] . may provide an alternative route to the same objective), particularly in the 
present era of precision numerical relativity where linear perturbation theory estimates may be 
compared directly with non-linear exact results. A key ingredient in the recipe for computing 
the degree of QNM ringing is the so-called 'QNM Excitation Factor': B ln = A^ ut) / \2u^-) . 
In a forthcoming paper, we will show how the expansion method can be combined with WKB 
techniques to find dA^/du and to compute excitation factors in the large-L limit. We will 
also show that the singular structure of the retarded Green function near the null cone may be 
deduce from the large-L asymptotics of QNM wavefunctions and excitation factors [581 137] . 

Further extensions of this work are possible. Firstly, we hope it would be a simple task 
to use the expansion method to explore the examples in Sec. [6] (or others) in greater depth, 
or to study the Regge pole frequencies and wavefunctions for these spacetimes. Secondly, 
there remains the question of how the expansion method could be extended to treat the 
Schwarzschild-AdS case (of interest in the context of the AdS/CFT conjecture [13]), where 
'reflecting' boundary conditions at r* — ► +oo are imposed. Thirdly, how might the expansion 
method be applied to the Kerr (rotating black hole) spacetime? Here, an added complication is 
that the angular separation constant depends on the frequency of the perturbation. However, 
it is known that a one-parameter (Q) family of unstable photon orbits at fixed (Boyer- 
Lindquist coordinate) r = r c (Q) exist [66J, so a geometrically-motivated line of enquiry may 
pay dividends. So far we have had some limited success in finding the large-Z asymptotics of 
the 'polar' (m = 0) mode of the scalar field. We hope to report on further progress in this 
direction in the near future. 
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